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ABSTRACT 

The spectra of baryons containing two heavy quarks test the form of the QQ 
potential through the spin-averaged masses and hyperfine splittings. The mass 
splittings in these spectra are calculated in a nonrelativistic potential model and 
the effects of varying the potential studied. The simple description in terms of 
light quark and pointlike diquark is not yet valid for realistic heavy quark masses. 



It is well-known that the quarkonium and T) spectra test the QCD potential 
between quark and antiquark. In this phenomenological flavor-independent 

potential and a QCD-derived hyperfine interaction can be used to describe the masses 
and splittings of the cc and bb mesons U|. Furthermore, such models allow predictions 
to be made for the unobserved states in the spectra |J and for the B c system ||. In the 
same way, the doubly-heavy baryons (those containing ccq or bbq) test the potential 
between two quarks. The one-gluon-exchange potential between heavy quarks differs 
from the quark- antiquark potential by a relative color factor 2. This relation does not 
hold at higher orders: the one-loop corrections to the QQ and QQ potentials are not 
equal ||]. The confining potential, which is determined phenomenologically, should 
be studied for QQ interactions as it has been for QQ. 

Tests of the QQ potential can be implemented not only in the spin-averaged 
spectra of the doubly-heavy baryons, but also by studying the hyperfine splittings of 
these particles. This is surprising, as the hyperfine interactions contain two types of 
terms: those describing the spin-spin interaction between the heavy quarks (~ l/wig), 
which depend on the QQ potential analogously to the quarkonium case, and those 
describing the spin-spin interactions between heavy and light quarks (~ l/mQm q ), 
which might be thought to depend on the light-heavy quark interaction and not 
significantly on the heavy-heavy quark potential. Nevertheless, the wave function of 
the light quark depends sufficiently on the heavy quark pair to provide information 
on the QQ binding. 

Baryon spectroscopy with non-relativistic potential models is relatively accurate, 
although not perfect. The method is clearly limited by its use of a non-relativistic 
treatment of the light quarks. Nonetheless, such models provide a relatively good 
description of the baryon spectra and stable predictions for the hyperfine splittings 
Sq — Eg and Eg — Ag for the baryons containing one heavy quark [0] . One potential 
choice which provides a good fit to the observed baryon spectrum is the power-law 
form ||: 

V(r tJ ) = W(A + B4 + -^-6(r ij )a i .a 1 ), (1) 
z i<j \ iihiiij j 

with 

j3 = 0.1, A = -8.337 GeV, B = 6.9923 GeV 1+/3 , C = 2.572 GeV 2 , 
m g = 0.3 GeV, m s = 0.6 GeV, m c = 1.895 GeV, m b = 5.255 GeV. 

The light and strange quark masses and the power (3 are held fixed at reasonable 
values, and the parameters A, B, and C are obtained from a fit to the masses of 
the N, A, and Q~. The charm and bottom quark masses are fit to the A c and 
This potential, when rescaled by the relative color factor of 2 between QQ and QQ 
one-gluon exchange, fits the J/\l/ and T spin-averaged spectra as well. Such power- 
law potentials are familiar from the description of QQ bound states and give the 
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choice of (3 used above. The value of C given above is however significantly larger 
than that obtained from a fit to the J/^/-r] c hyperfine splitting, which is 1.172 ||. 

Although this potential fits the known baryon spectrum, it is not convenient for 
study of the effects on the spectra of variation of the QQ potential. In particular, it is 
interesting to be able to vary separately the Coulomb-like small distance contribution 
to the potential and the long-distance confining term. The potential between light 
and heavy quarks is therefore left as in Eq. (p]), while the QQ interaction is replaced 
with the Cornell form [Tcfl : 



VQQ(r) = U-- + ar + c + -Z r 5(r)ff 1 -a 2 ), (2) 
z. \ r m Q j 

a = 0.2005 GeV 2 , k = 0.4274, c = -1.0455 GeV. 

The parameters a, k, and c are obtained from a fit to the J/ty and T spectra, with the 
quark masses and C held fixed to the values given above. The effects on the hyperfine 
splittings of small variations in a and k are then studied. These variations, as well as 
the power-law type potential given above, are plotted in Figure 1 for comparison. 

In the limit of heavy quark mass, the heavy QQ pair is more easily excited to 
higher-energy states than the light quark. In the harmonic oscillator potential, for 
example, the ratio of the excitation energies for light and heavy degrees of freedom is 
\J (2mh + mi)/3rri[ [|TTJ. The states considered here are those containing two identical 



heavy quarks (cc and bb rather than be) in a relative S-wave state, so that Fermi 
statistics require that the heavy quarks be in a spin-1 (symmetric) state. In the 
standard notation the lowest-energy states for QQu or QQd are Hqq (^qq) f° r s Pi n 
1/2 (3/2), and for QQs they are denoted by Qqq (^qq )• The heavy- heavy hyperfine 
terms cancel in the difference between spin-3/2 and 1/2 states as the expectation 
value of the spin-spin operator, < B\ • a 2 >, is identical for these two cases. The 
hyperfine splittings discussed below thus originate solely from the light-heavy quark 
interactions. Three different mass splittings are considered here for each combination 
of heavy and light quark flavors. The first is the difference between the centers-of- 
mass of the states in which the QQ pair is in 2S and IS states. The other two are 
the hyperfine splittings between the spin-3/2 and spin- 1/2 baryons for the 2S and IS 
QQ configurations. 

Once the potentials are set, the solution of the three-body problem can be ap- 
proached using a variety of numerical methods [f[T]| . The Born-Oppenheimer method 
is an efficient solution method for the case of two heavy and one light quark. The 
wavefunction is split into heavy- and light-quark degrees of freedom, 

*(p,A)=5>n(p)/n(p,A), (3) 

n 

where p is the distance between the two heavy quarks (mi and mz), and A = 



3 



(r*3 — Ri2)\Jrn 3 (mi + m2)/ /J,i2(m 1 +m 2 +m 3 ) is proportional to the distance between 

the light quark m 3 and the center-of-mass R\ 2 of the heavy-quark pair. This choice of 
variables separates the non-relativistic kinetic energy terms: p\/(2m\) + p\j '(2m 2 ) + 
p\/{2m 3 ) = \p 2 p + p 2 x }/ (2fi 12 ), where /ii 2 is the reduced mass of the two heavy quarks. 
The dependence on the light-quark mass has been absorbed into the definition of A. 
The light quark wavefunction / is found for fixed p 



2fi 



1 Pl + V(r 13 ) + V(r 23 ) 
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fn{p,X) = e n (p)f n (p,X), (4) 



where is the distance between the light quark and heavy quark i. The energy of 
the QQq state is approximated by: 



' vl + Vqq(p) + e (p)+ < fo\^—pl\fo > 
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Up) = eMp)- (5) 



Keeping only the zeroth term in the expansion of \I/ corresponds to the so-called 
"uncoupled adiabatic" [TI| approximation. 



The Born-Oppenheimer method integrates over the distances between the heavy 
quarks in the second step (Eq. (|)), so that the above approximation differs from 
that in which the heavy-quark pair behaves as a spin-one pointlike diquark, as it 
would for infinitely heavy quark mass. Because this method does not depend on 
the limit of infinite tuq, it works well already for the charmed baryons. In the limit 
7^1,2 m 3 these equations approach the expected form: The distances r i3 in Eq. (||) 
are independent of the quark masses, and the term p\j2\iyi — ► p\/2m 3) so that the 
light-quark wavefunction is independent of the distance p between the heavy quarks. 
Eq. (Bj) then reduces to the usual non-relativistic two-body equation. In the limit of 
infinite tuq, the QQq spectra can be simply related to the meson spectra 0. The 
splittings obtained in this way are however significantly smaller than those obtained 
from the potential model calculation, and the masses of the spin-averaged states 
somewhat higher. This limit gives the values: M=* c = 3.811 GeV, M= cc = 3.742 GeV, 
M 3 * bb = 10.484 GeV, and M Ebb = 10.460 GeV §. 

The Coulomb and linear terms in the QQ potential are now varied by introducing 
parameters x a and Xk'- 

1 / k C \ 

v QQ( r ) = - \-x k - + x a ar + c-\ T 5(r)a x ■ a 2 ), (6) 

z \ r m Q J 

which are varied between 0.8 and 1.2 to determine the effect on the mass differences 
M(2S)-M(1S), AM(1S), and AM(2S). Variations in c do not affect the splittings, only 
the masses themselves. Variations in C obviously cause proportional changes in the 
hyperfine splittings. The effect of C on the 2S-1S mass difference comes from the 
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1/mg hyperfine terms only, and is small. These results of the x a and Xk variations 
are presented in Table [l], along with the percent change 100(AM(sj = 1.2) — AM(xj = 
0.8))/AM(ii = 1). The values for x a = = 1 are in agreement with previous 
calculations |L4| when the power-law type potential is used for the QQ interaction. 

The hyperfine splittings are analogous to those in a heavy-light meson system, 
so that in the limit of infinitely heavy itlq they should approach S/tuq, where 5 is 
independent of heavy-quark flavor. The mass splittings in the table give 5(ccq,lS) = 
0.242 and <5(bbq,lS) = 0.305. These states are not yet heavy enough to be accurately 
described by the simple heavy-quark limit. The values of 5 extracted for these and 
for heavier masses are shown in Fig. 2, where 5 may be seen to approach the heavy- 
mass limit for mg> 15-20 GeV. This is in contrast to what is known from meson 
spectroscopy, where the simple l/mg description of the splittings works even for the 
strange quark. Although the masses of the diquarks are quite large (~ 2mg), the radii 
for these states are also large for realistic quark masses, so that the simple application 
of the infinite-mass limit is not yet valid here. Typical radii < r > for QQ states are 
2.4 GeV -1 for cc IS, and 1.5 GeV -1 for bb IS. For QQ states, typical radii are 1.7 
Ge^ 1 for cc, and 1.0 GeV -1 for bb. 

The percent changes in the mass splittings are given in Fig. 3 for QQq. The 
sensitivity of the QQs states to changes in the potential is almost exactly that of 
the QQq states. M(2S)-M(1S) is more sensitive to changes in the QQ potential 
than either of the hyperfine splittings. This mass difference is essentially the mass 
splitting between the IS and 2S states of the two-body QQ potential, so that the 
strong dependence here is unsurprising. The dependence of M(2S)-M(1S) on x fc is 
about one-quarter(half) that on x a for the cc (bb) states. The bb pair is more closely 
bound than the cc so that these states depend more strongly on the Coulomb (xk) 
part of the potential. 

The hyperfine splitting in the IS states, AM(1S), depends more strongly on x a 
than on x^ for the cc case, and on x^ than on x a for bb. The strong dependence of 
the bb IS splitting on the Coulomb term in the potential reflects the closer binding 
of the b quarks. On the other hand, AM(2S) depends in both cases more strongly on 
x a than on x^. This is a reflection of the larger radius (< r >~3.0 GeV -1 ) of the bb 
2S state. 

In conclusion, the spectra of doubly-heavy baryons have been studied in a non- 
relativistic potential model, with emphasis on the sensitivity of the mass splittings of 
these states to variations in the interaction between the two heavy quarks. Although 
the hyperfine splittings are between the two possible spin states of the light quark 
relative to a heavy spin-1 system, it was found that dependence on the QQ interac- 
tion remains. This dependence is however rather small. The 2S-1S mass difference 
is more sensitive to variations in the potential than the hyperfine splittings. A com- 
parison of the splittings between the center-of-mass of the 2S and IS states and the 
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hyperfine splittings in these levels exhibits different dependences on the strengths of 
the separate parts of the potential (Coulomb vs. confining terms). More importantly, 
it is demonstrated that the use of the infinitely heavy mg limit, where the heavy 
quarks form a pointlike diquark and the baryon as a whole a meson-like system, is 
not yet valid for realistic heavy quark masses. 
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Figures 

Fig. 1. The QQ potentials used in this paper. Solid lines denote the power-law 
(Eq. (0)) and Coulomb-plus-linear potentials (Eq. (|2|)) without any varia- 
tion. The dot-dashed (dotted) lines show the effect of the variation of Xk 
(x a ) between 0.8 and 1.2. 

Fig. 2. The quantity 5 = mgAM(lS) as a function of heavy-quark mass. 

Fig. 3 . The percent changes in the mass splittings 8A4 (x j) = 1 00 (.M (xj) /.M (1) — 1 ) • 
The labels 1, 2, 3 refer to M(2S)-M(1S), AM(1S), and AM(2S), respectively, 
while a and k indicate variations in x a or x^. 







M(2S)-M(1S) 


AM(1S) 


AM(2S) 


ccq 


0.80 


0.394 


0.125 


0.0802 




1.00 


0.431 


0.128 


0.0832 




1.20 


0.467 


0.130 


0.0861 




% change 


16.9 


4.2 


7.0 


bbq 


0.80 


0.280 


0.0571 


0.0417 




1.00 


0.304 


0.0581 


0.0430 




1.20 


0.328 


0.0590 


0.0441 




% change 


15.7 


3.3 


5.6 












ccq 


0.80 


0.423 


0.126 


0.0828 




1.00 


0.431 


0.128 


0.0832 




1.20 


0.440 


0.129 


0.0837 




% change 


3.8 


2.3 


1.1 


bbq 


0.80 


0.294 


0.0565 


0.0425 




1.00 


0.304 


0.0581 


0.0430 




1.20 


0.316 


0.0593 


0.0435 




% change 


7.4 


4.8 


2.4 



Table 1: Mass splittings (GeV) in the ccq and bbq systems as a function of the 
parameters x a and xy.. 
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